Kis transition form factors 



Chueng-Ryong Ji and Pieter Maris 
Department of Physics, North Carolina State University, Raleigh, NC 27695-8202 

(February 6, 2008) 

The rainbow truncation of the quark Dyson-Schwinger equation is combined with the ladder 
Bethe-Salpeter equation for the meson bound state amplitudes and the dressed quark-W vertex in 
a manifestly covariant calculation of the Ki^ transition form factors and decay width in impulse 
approximation. With model gluon parameters previously fixed by the chiral condensate, the pion 
mass and decay constant, and the kaon mass, our results for the Ki-j, form factors and the kaon 
semileptonic decay width are in good agreement with the experimental data. 
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^ . I. INTRODUCTION 

D ' The central unknown quantity required for reliable calculations of weak decay amplitudes are the hadronic matrix el- 
pLn \ ements. In this respect, the semileptonic Ki^^ decay is an interesting process: pseudoscalar mesons are well-understood 

• as Goldstone bosons, only the vector part of the weak interaction is involved, and there are experimental data available 
^SJ , which allows one to judge various theoretical (model) calculations. These calculations generally fall into two different 

classes: effective theories using meson degrees of freedom, e.g. chiral perturbation theory 0J^, and models employing 
' quark and gluon degrees of freedom, e.g. constituent quark models |3|-||] and the Dyson-Schwinger approach 1^-^. 
^ , In particular, the light-front approach has been a popular framework for a Hamiltonian approach to analyze exclusive 
hadronic processes , and the pseudoscalar-to-pseudoscalar weak transition form factors have been studied within 
a light-front constituent quark model in Ref. iQ,^. One complication of the light-front formalism is that in the 
analysis of timelike exclusive processes one needs to take into account light-front nonvalence contributions, which 
I are absent in a manifestly Poincare invariant approach. Recently, an effective treatment has been presented [^0| to 
T-H . incorporate such contributions, and a systematic program has been laid out lITl]] to take into account the nonvalence 
■ contributions. However, the systematic program explicitly requires all higher Fock-state components while there has 
: been relatively little progress in computing the basic wave- functions of hadrons from first principles. Furthermore, the 

• (nonperturbative) dressing of the quark-W vertex and the dynamical chiral symmetry breaking are yet underdeveloped 
i-H \ aspects in this approach. In particular the latter aspect needs to be seriously considered for processes involving pions 

^ ' and kaons, which are the pseudo-Goldstone bosons associated with dynamical chiral symmetry breaking p^ , p^ . 
(— ( In this work we use the set of Dyson-Schwinger equations [DSEs] to calculate the Ki^i transition form factors. 

• • . The DSEs provide us with a manifestly covariant approach which is consistent with dynamical chiral symmetry 
breaking jl^ , electromagnetic current conservation Jl5[ , and quark and gluon confinement ]l6t . For recent reviews on 



in 
o 



, the DSE and its application to hadron physics, see Rcfs. |17 1^. Our calculation of the Ki^ transition form factors is 
^ ' analogous to recent calculations of the pion and kaon electromagnetic form factors in impulse approximation | [l5[ | . Here 
we demonstrate that in the same framework, with a consistent dressing of the quark-W vertex, we can also describe 
the semileptonic decay of both the neutral and charged kaons, without any readjustment of the model parameters. In 
Sec. II we review the formulation that underlies a description of the meson form factors within a modeling of QCD 
through the DSEs, and discuss the details of the model. In Sec. Ill we present our numerical results for the form 
factors and decay widths. A discussion of our results is given in Sec. IV. 

II. MESON FORM FACTORS WITHIN THE DSE APPROACH 

The matrix elements describing kaon semileptonic decays are 

jf {P,Q) = {^\p)\s^^u\K+{k)) 

= ^ {f+{~Q^)P^ + /-(-Q')Qp) , (1) 

for the charged kaon, and for the neutral kaon 
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{P,Q) = {7:-{p)\-sj,u\K'{k)) 

= f+{-Q')P,, + f-{-Q')Q^., (2) 

where = (p + fc)^ and = (fc — p)^, with P ■ Q = ni^ — for on-shcU pions and kaons[]. Ahernatively, we can 
decompose into its transverse and longitudinal components 

jf{P.Q) = f+i-Q^) T^u{Q)Pu + /o(-O') ^ Q,. , (3) 

with T^u the transverse projection operator 

T,.(Q)=(V-^) • (4) 

In the isospin-symmetric limit, which we employ here, the form factors f±{t = — Q^) and /o(i) are the same for the 
K+ and the 

In impulse approximation, these matrix elements are given by 

jf {P, Q)=N,j ^ 'Tr[S\q) VjUq, q - k)S^{q - k) iVf'^iq -~k,q-p) S^{q - p) ff{q - p, q)] , (5) 

where S-^ is the dressed quark propagator with flavor index /, T^(^j(-j the pion (kaon) Bethe-Salpeter amplitude [BSA], 
and rj^"^ the dressed suW vertex, each satisfying their own DSE. Note that the coupling constant and the CKM 
matrix element Vus are removed from the definition of the quark-W vertex. 

A. Dyson— Schwinger Equations 

The DSE for the rcnormalizcd quark propagator in Euclidean space is 

Sip)-'=zZ2^ + Z^mi^i) + Z, p-0-^g^D^,{k)^^^S{q)Tl{q,p), (6) 

where Dp_^{k) is the dressed-gluon propagator, T1{q]p) the dresscd-quark-giuon vertex, and k = p — q. The most 
general solution of Eq. (||) has the form S{p)~^ = i^A{p'^) + B{p'^) and is renormahzed at spacelike jj? according to 
A{iJ?) = 1 and B{iJ?) = m{n) with m(/z) the current quark mass. The notation represents a translationally invari- 
ant regularization of the integral with A the regularization mass-scale; at the end of all calculations this regularization 
scale can be removed. 

Mesons can be studied by solving the homogeneous Bethe-Salpeter equation [BSE] for bound states, with a 
and b flavor indices, 

r^^(p+,p_)= p-^^K{p,q;Q)®S-{q+)T'^^{q+,q^)S\q^), (7) 

where — P + vQ ^^.d p_ = p — {1 — ri)Q are the outgoing and incoming quark momenta respectively, and similarly 
for q±. The kernel K is the renormahzed, amputated qq scattering kernel that is irreducible with respect to a pair of 
qq lines. This equation has solutions at discrete values of = —mfj, where niH is the meson mass. Together with 
the canonical normalization condition for qq bound states, it completely determines Th, the bound state BSA. The 
different types of mesons, such as (pseudo-)scalar, (axial-)vector, and tensor mesons, are characterized by different 
Dirac structures. The most general decomposition for pseudoscalar bound states is ||l3| 

Tps{k+,k^) = j5[iE{P;k ■ Q;r,) + Q Fik^;k ■ Q;tj) + ^G{k^-k ■ Q-r^) + a^.Q^k, H{k^-k ■ Q-r^)] , (8) 



'^We employ a Euclidean space formulation with < for timelike vectors, 7^,71/ = 25^^, and 7^ = 7m- 



2 



where the invariant ampHtudes E, F, G and H are Lorentz scalar functions of fc^ and k-Q. Note that these amphtudes 
expHcitly depend on the momentum partitioning parameter 77. However, so long as Poincare invariance is respected, 
the resulting physical observables are independent of this parameter jTs]] . 

In order to describe weak (and electromagnetic) form factors one also needs the dressed qqW (and qq^) vertices. 
These vertices satisfy an inhomogeneous BSE: e.g. the suW vertex rjj*^(p+,p_) satisfies 

Tf^ip+.p.) = (7p - 7m75) + K(p, q; Q) ® S^{q+) rf'^(g+, q-) S^{q-) . (9) 



B. Model Truncation 



To solve the BSE, we use a ladder truncation, with an effective quark-antiquark interaction that reduces to the 
perturbative running coupling at large momenta p^ , p^ 

K{p, q; P) ^ -g{e) i?fr-(fc)^7^ ® ^7, , (10) 

where Z)|j^''(fc — p ~ q) is the free gluon propagator in Landau gauge. The corresponding rainbow truncation of the 
quark DSE, Eq. (1), is given by ri{q,p) 7^AV2 together with Zig'^D^,^{k) g{k^)D^^;;''{k). This combination of 
rainbow and ladder truncation preserves both the vector and axial- vector Ward-Takahashi identities. This ensures that 
pions are (almost) massless Goldstone bosons associated with dynamical chiral symmetry breaking p3[ ; in combination 
with impulse approximation for meson form factors, Eq. (^, it also ensures current conservation. Furthermore, this 
truncation was found to be particularly suitable for the flavor octet pseudoscalar and vector mesons since the next- 
order contributions in a quark-giuon skeleton graph expansion, have a significant amount of cancellation between 
repulsive and attractive corrections [ pO| . 

The model is completely specified once a form is chosen for the "effective coupling" Q{k^)- We employ the Ansatz |l^] 
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(11) 



with 7,„ — 12/(33 — ^Nf) and T{s) = (1 — exp j^)/s. The ultraviolet behavior is chosen to be that of the QCD 

running coupling a(fc^); the ladder-rainbow truncation then generates the correct perturbative QCD structure of 
the DSE-BSE system of equations. The first term implements the strong infrared enhancement in the region 

< A:^ < 1 GeV^ phenomenologically required [|T] to produce a realistic value for the chiral condensate. We use 
mt = 0.5 GeV, r = e^ — 1, Nf — 4, Aqcd = 0.234 GeV, and a renormalization scale /i = 19GcV which is well into 
the perturbative domain [^3|jl^. The remaining parameters, uj = 0.4 GeV and D = 0.93 GeV^ along with the quark 
masses, are fitted to give a good description of the chiral condensate, m^/if and /jr. 

Within this model, the quark propagator reduces to the perturbative propagator in the ultraviolet region. However, 
in the infrared region both the wave function renormalization Zij?^ = XjAijP'^ and the dynamical mass function 
MijP') — B{p'^)/A{p'^) deviate significantly from the perturbative behavior, due to chiral symmetry breaking. It is 
interesting to note that the typical results obtained using the rainbow DSE, a significant enhancement of M(p^) below 

1 GeV^ and also an enhancement of A{p'^), have recently been confirmed by lattice simulations |22); with the present 
ladder DSE model, the functions Z and M are in semiquantitative agreement with the forms obtained in lattice 
simulations ||2^ . 

The vector meson masses and electroweak decay constants one obtains in this model are in good agreement with 
experiments |l9). Without any readjustment of the parameters, one gets remarkable agreement with the most recent 
Jlab data for F^^ ||2^ , and also for other electromagnetic charge radii and form factors [|l5|,^ . The strong decays of 
the vector mesons into a pair of pseudoscalar mesons are also well-described within this model |p6|| . Here we apply 
the same approach to the kaon semileptonic decay. 



III. NUMERICAL RESULTS 



To calculate the form factors f±{t), we start by solving the quark DSE in rainbow approximation, and subsequently 
use its solution to solve the homogeneous BSE, Eq. (0), for the pseudoscalar bound states, and the inhomogeneous 
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BSE, Eq. (||), for the quark-W vertex, using the ladder truncation, Eq. (|o|). With these BS amphtudes and quark 
propagators, we then calculate the form factors in impulse approximation using Eq. (I). In the SU(3) flavor Hmit, 
fj^ — F-^, the pion electromagnetic form factor, and /_ ~ 0. We therefore expect /+ to be of order one for small t. 



On the other hand, /_, which is a measure of the constituent mass ratio Ms/M^ 
smaller. This expectation is indeed confirmed by our calculations, see Fig. ^. 
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FIG. 1. Our results for (solid dots) with linear fit (dashed line), and f^{t) (solid diamonds) with a constant fit 

(dot-dashed hue) in the physical region. 

In the physical region, mf < t < (mK^m^)'^ ~ 0.13 GeV^, our results for f±{t) are almost linear. The experimental 
data are often analyzed in terms of /+ and /o, using a linear parametrisation 



/(0 = /(0)(l + 



(12) 



Such a linear parametrisation for both /+ and /o implies that /_ is independent of t, since /o is related to f± via 

t 



m - U{t) + 



(13) 



To within a few percent, /_ is indeed almost constant in the physical region, < t < 0.13 GeV^, as can be seen from 
Fig. ^. Our results are summarized in Table |, where we have included a slope parameter A for /_ {t) as well. 

Our results agree quite well with the available data for f+{t), see Fig. |. Also the results for other observables 
compare reasonably well with the experimental data, given the error bars, and the differences between the form factor 
parameters extracted from the charged and the neutral kaon decay, see Table ||. The partial decay width for K n eve 
and K n /it/^ is obtained by integrating the decay rate 



dr 



Gf \Vus\^ 

247r3 



where 



(14) 



K{t) 



2mK 



{{ml + ml - - imlml) 



2\2 



and mi is the lepton mass. Again, we find good agreement with the data, see Table |[ 



(15) 
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FIG. 2. Our results for f+{t), normalized to ,/+(0) = 1, compared to the experimental results 



For comparison, we have also included the results from some other model calculations ||7|,|10| and from 1-loop chiral 
perturbation theory Our results compare quite well with chiral perturbation theory, not only for the slope 

parameters, but also for the curvature of the form factors. A prediction from current algebra is the value of /o(i) at 



the Callan-Treiman point, t 



'if ■ 



-mi, namely /if/A = 1.22 



The reason we get foimj^ — mi) = 1.18 instead 



of 1.22 is related to our value of fx'- we obtain Jk = 109 MeV in our model |1^J19|], compared to the experimental 
value 113 MeV. Because our approach satisfies constraints coming from current conservation and the axial- vector 
Ward-Takahashi identity, it is not surprising that we indeed agree with these predictions from chiral perturbation 
theory. 
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FIG. 3. Our results for f+{t) (solid dots) with monopole fit (solid curve), and fo{t) (open squares) with a constant plus 
monopole fit (dashed curve). 



Outside the physical region, f±{t) and fo{t) are clearly not linear, as shown in Fig. C 
our result for f+{t) is very well approximated by a monopole, with monopole mass m 



. In the region we considered, 
= 0.87 GeV^. The reason for 

this monopole behavior is easy to understand: the vertex function p^"^ has a resonance peak in its transverse vector 



components at the K* mass, because the homogeneous version of Eq. (O) has a solution at — 



'if 



, corresponding 



to the K* vector meson. The longitudinal part of the quark-W vertex has a resonance peak at the location of the 
us scalar bound state. We tentatively identify this with the Kq(1A3Q) meson; however, in contrast to vector states, 
the use of ladder approximation is expected to be inadequate for scalars [Q. Near these bound states, the vertex 
behaves like 



W 



(p+,p_;g)^r-g"i-(p+,p„;Q) 



'if* 



r^«(p+,P-;Q) 



(16) 



where F-'^ and F^o are the bound state BSAs, and ry and rs the residues of the vertex at these poles In our 
model we have actual poles in the quark-W vertex, rather than resonance peaks, because we have not included meson 
loop corrections in our kernel, and therefore we do not generate a width for the K* nor the K^. In the spacelike 
region and at small timelike t — —Q^ we expect these corrections to be small. However, they will modify our results 
significantly close to the resonance peak. 
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As can be seen from Eq. these poles in F*"^ lead a similar behavior of near these poles. From Eq. (^) it 
is clear that the vector meson pole leads to a pole in /+(t), but not in fo{t), whereas the scalar meson induces a pole 
in fo{t) but not in f+{t) 

foi-Q') = , (17) 

Note that f-{t) will in general be sensitive to both the scalar and the vector meson bound state. 

The present model has a vector bound state at t = 0.876 GeV^ (cf. the experimental mass ra\, = 0.796 GeV^) pi. 
Thus the monopole behavior of the form factor f+{t) in Fig. ^ simply reflects the existence of this pole. A similar 
behavior was found in the pion and kaon electromagnetic form factors [T5| ] . The physical /+ {t) can be expected to rise 
more slowly with t than our present calculations due to meson loop corrections, and will develop a nonzero imaginary 
part above the threshold for intermediate tt K states. 

The scalar form factor /o(i) does also exhibit the expected pole behavior, and we are able to identify the scalar and 
the vector poles separately, see Fig. ^. The lowest us scalar bound state in this model has a mass of m? = 0.80 GeV^, 
which is rather low compared to the i4rQ(1430) with a mass m? = 2.04 GeV^; on the other hand, it could be an 
indication that there exists a light scalar resonance k(900), as has been speculated 133]. Something similar happens in 
in the u/d quark sector: in rainbow-ladder truncation, one typically finds a scalar ujd bound state around m = 0.6 to 
0.7 GeV [^-^, just below the p mass. This bound state could either correspond to a broad ct, or to the /o and/or aq, 
in which case the calculated mass is about 30 to 40% too low. It is known that the leading perturbative corrections to 
the ladder kernel cancel to a large extent in both the pseudoscalar and vector channel, but not in the scalar channel; 
we therefore expect in the scalar channel significant repulsive corrections to the ladder kernel which could increase 
the mass of the scalar bound state We also expect meson loop corrections to be more important in the scalar 
channel than in the vector channel. Therefore, our calculation for fo{t) should not be trusted quantitatively in the 
timelike region beyond the Callan-Treiman point, t = 0.23 GeV^. 

In the spacelike asymptotic region /+ and /o seem to behave differently: a monopole fit does not work for /q. In 
order to fit /o, we either need to add a nonzero constant, or use at least two monopoles. Simple power counting 
indicates that jj^ scales like l/Q; inserting this behavior in Eqs. ( p^ ) and ( p^ implies that f+{t) falls like but 
/o(t) goes to a constant (up to logarithmic corrections) at large spacelike t (remember that P scales like Q in the 
asymptotic region). This behavior of /+ agrees with the pQCD predictions for the electromagnetic pion and kaon 



form factor |36|; the behavior of /o can be understood if one realizes that the combination /o(t) ""-^ '"'^ does fall like 



1/^ if /o(i) goes to a constant. 



IV. DISCUSSION 



Using the rainbow-ladder truncation of the set of DSEs with a model for the effective quark-antiquark interaction 
that has been fitted to the chiral condensate and /tt, we study the Ki^ decay in impulse approximation. Our results, 
both for the form factors and for the decay width, are in good agreement with experimental data, without any 
readjustment of the parameters; they also compare quite well with chiral perturbation theory. Note however that 
in chiral perturbation theory the pion charge radius is used as input, in order to fix a low-energy constant, which is 
important for these form factors as well, in particular for A_|_. In our calculation, the only model parameters are in 
the infrared behavior of the effective quark-antiquark interaction, which were fitted to the chiral condensate and the 
pion decay constant [Q; the pion charge radius follows from a calculation similar to the one presented here for the 
Ka decay 

Our approach is based on previous work by Kalinovsky et al. 0, and the results are quite similar. However, an 
important difference with this earlier study is that we dress the quark-W vertex by solving the inhomogeneous BSE 
for this vertex. The main advantage of doing so is that we thus automatically include effects coming from intermediate 
vector and scalar mesons. Another difference is that here we use actual solutions of a DSE for the propagators and 
BS amplitudes in Eq. (|^), rather than phcnomcnological parametrisations, which reduces the number of parameters 
in the calculation. 

We have demonstrated explicitly that both /+(t) and /o(t) exhibit resonance peaks in the timelike region due to 
the existence of vector and scalar bound states respectively. The effects of meson loops are not included in our present 
calculation, which is why we find a pole behavior rather than a resonance peak. In Ref. [Fj it was already demonstrated 
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that non-analytic effects from such loops could contribute significantly to the behavior of the form factors, not only 
beyond the threshold for ttK production, but also in the physical region. We hope to be able to incorporate these 
effects in future work. 

Our results are also similar to those obtained recently by Ji and Choi in a light-front calculation, at least 
for f+(t), which dominates this decay, even though details of the calculation are quite different. In light-front 
calculations of timclike processes one has to include particle-number-nonconserving Fock states to recover the Lorentz 
covariance, which complicates the calculation Since our approach is manifestly covariant, such contributions 

are automatically included in our impulse approximation. Another difference is that we use momentum-dependent 
quark self-energies, consistent with dynamical chiral symmetry breaking, whereas in Ref. [|lO| constituent quarks with 
fixed masses are used. Recently, the effects of a running quark mass (instead of a fixed constituent mass) have been 
explored in a light-front calculation of the pion form factor |37| . It would be interesting to see the effect of such a 
running mass on the K13 form factors, in particular on f~{t), which in general appears to be quite sensitive to details 
of the calculation. 
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our 


experiment 


xPT iii 


other theory 


light cone @ 




calc 


K+ 


K" 


DSE model § 


/+(0) 


0.964 






0.977 


0.98 


0.962 


A+(e3) 


0.027 


.0276 ± .0021 


.0288 ± .0015 


0.031 


0.028 


0.026 


A+(^3) 


0.027 


.031 ± .008 


.034 ± .005 


0.031 


0.029 


0.026 


-/-(o) 


0.10 






0.16 


0.24 




A- 


0.03 








0.023 




-e = -/-(o)//+(o) 


0.11 


0.31 ± .15 


0.11 ± .09 


0.17 


0.25 


0.01 


Ao 


0.018 


.006 ± .007 


.025 ± .006 


0.017 


0.007 


0.025 


M-rriK - ml) 


1.18 






1.22 


1.18 






7.38 


3.89 


7.50 






7.3 


r(^M3) 


4.90 


2.57 


5.26 






4.92 



TABLE I. Our results, compared to data ^g] where available and some other calculations. The partial decay width is in 
10® s~^; because — /\/2 in the exact isospin limit, the experimental value for the K'^ partial decay width should be 
multiplied by a factor of two, in order to compare it with our calculation and with the K*^ partial decay width. 
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